Magnetic reversals in a simple model of MHD 
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We study a simple magnetohydrodynamical approach in which hydrodynamics and MHD turbu- 
lence are coupled in a shell model, with given dynamo constrains in the large scales. We consider 
the case of a low Prandtl number fluid for which the inertial range of the velocity fleld is much wider 
than that of the magnetic field. Random reversals of the magnetic field are observed and it shown 
that the magnetic field has a non trivial evolution - linked to the nature of the hydrodynamics 
turbulence. 
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Observations show that natural dynamos are intrinsi- 
cally dynamical. Complex magnetic field evolutions have 
been reported for many systems, including the Sun and 
the Earth [l[. Formally, the coupled set of momentum 
and induction equations are invariant under the trans- 
form: (u,B) (u, — B) so that states with opposite po- 
larities can be generated from the same velocity field (u 
and B are respectively the velocity and magnetic fields). 
In the case of the geodynamo, polarity switches are called 
reversals [l[ and occur at very irregular time intervals ■ 
Such reversals have been observed recently in labora- 
tory experiments using liquid metals, in arrangements 
where the dynamo cycle is either favored artificially Q 
or stems entirely from the fiuid motions 0, Q . In these 
laboratory experiments, as also presumably in the Earth 
core, the ratio of the magnetic diffusivity to the viscos- 
ity of the fluid (magnetic Prandtl number Pa/) is quite 
small. As a result, the kinetic Reynolds number Rv of 
the flow is very high because its magnetic Reynolds num- 
ber Rm = RvPm needs to be large enough so that the 
stretching of magnetic fields lines balances the Joule dis- 
sipation. Hence, the dynamo process develops over a 
turbulent background and in this context, it is often con- 
sidered as a problem of 'bifurcation in the presence of 
noise'. For the dynamo instability, the effect of noise en- 
ters both additively and multiplicatively, a situation for 
which a complete theory is not currently available. Some 
specific features have been ascribed to its onset (e.g. bi- 
furcation via an on-off scenario Q) and to its dynam- 
ics . Turbulence also implies that processes occur over 
an extended range of scales; however, in a low magnetic 
Prandtl number fluid the hydrodynamic range of scales is 
much wider than the magnetic one. In laboratory exper- 
iments, the induction processes that participate in the 
dynamo cycle involve the action of large scale velocity 
gradients [3, d, Q , with also possible contributions of ve- 
locity fluctuations at small scales 1^, 11. 12|. 



bulent fluctuations (as opposed to 'noise') in a dynamo 
instabilit y. The models stems from the approach intro- 
duced in[l3| for the hydrodynamic studies. Magnetic 
field reversals are observed above onset and we detail 
their characteristics. 




Building upon the above observations, we propose here 
a simple model which incorporates hydromagnetic tur- 



10 100 1000 10000 100000 1e+06 



FIG. 1: Main figure: the behavior of (I-B2P) (red triangles) 
anddBaP) ( green circles) as a function of Um for fixed value of 
V = 10~^. The blue line with dot corresponds to the solutions 
of ([5]) . Upper insert: energy spectra for _B„ (green triangles) 
and Un (red circles) corresponding to the = 0.001. 

Lower insert: the amount of magnetic dissipation (red tri- 
angles) tB = (E„z^„ifc^|_B„p) and the dissipation due to the 
large scale term e™ = am(|_B2p(S|r — BIJ) (green triangles). 



We consider an 'energy cascade' model i.e. a shell 
model aimed at reproducing few of the relevant charac- 
teristic features of the statistical properties of the Navier- 
Stokes equations [31 ■ In a shell models, the basic vari- 
ables describing the 'velocity field' at scale r„ = 2~"ro = 
k~^, is a complex number zi„ satisfying a suitable set of 
non linear equations. There are many version of shell 
models which have been introduced in literature. Here 
we choose the one referred to as Sabra shell model. Let us 
remark that the statistical properties of intermittent fiuc- 
tuations, computed either using shell variables or the in- 
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stantaneous rate of energy dissipation, are in close quali- 
tative and quantitative agreement with those measured in 
laboratory experiments, for homogeneous and isotropic 
turbulence IJ] . MHD shell model ~ introduced in [l5| - 



allow a description of turbulence at low magnetic Prandtl 
number since the steps of both cascades can be freely ad- 
justed 3, 17|. Although geometrical features are lost, 
this is a clear advantage over 3D simulations [3, fiot - 
We consider here a formulation extended from the Sabra 
hydrodynamic shell model: 



dun 
"dt 
dBn 
dt 

where 



-($„(U, U) - $„(B, B)) - Vklun + fn , (1) 

^($„(u, B) - $„(B, u)) - u,nklB„ , (2) 



$„(u, w) = fc„+i[(l + (5)u„+2W*+i + (2 - (5X+iu;„+2] 

+ fc„[(l - 2S)ul^j^Wn+l - (1 + (5)lt„+lW*_J 

+ fc„_i[(2 - (5)w„_iU'„_2 + (1 - 2S)Un-2Wn-l] , (3) 

for which following we chose 6 ~ —0.4. For this 
value of 6, the S'abra model is known to show statis- 
tical properties (i.e. anomalous scaling) close to the 
ones observed in homogenous and isotropic turbulence. 
The model, without forcing and dissipation, conserve 
the kinetic energy S„|m„P, the magnetic energy Eji|B„p 
and the hclicity i?e(I]„u„i?*). In the same limit, the 
model has a U{1) symmetry corresponding to a phase 
change exjp{i6) in both complex variables u„ and i?„. 
The quantity $„(?j,w) is the shell model version of the 
transport term ij\/w. The forcing term /„ is given by 
/„ = Sinfo/ul, i.e. we force with a constant power in- 
jection in the large scale. We want to introduce in eq. 
([2]) an extra (large scale) term aimed at producing two 
statistically stationary equilibrium solutions for the mag- 
netic field. For this purpose, we add to the r.h.s. of ^ 
an extra term M2{B2), namely for n ~ 2 eq.® becomes: 
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(4) 



where F2{u,B) is a short hand notation for 
i/3($2(it,5) - $2(5,"))- The term M2{B2) is 
chosen with two requirements: 1) it must break the 
U{1) symmetry; 2) it must introduce a large scale dis- 
sipation needed to equilibrate the large scale magnetic 
field production. There are many possible ways to 
satisfy these two requirements. Here we simply choose 
M2{B2) = flm-Bi- We argue, see the discussion at 
the end of this letter, that the two requirements are a 
necessary condition to observe large scale equilibration. 
From a physical point of view, symmetry breaking 
also occurs in real dynamos since the magnetic field is 
directed in one preferential direction which changes sign 
during a reversal. Thus symmetry breaking is a generic 
feature which wc introduce in our model by prescribing 



some large scale geometrical constrain. On the other 
hand, large scale dissipation must be responsible of the 
equilibration mechanism of the large scale field. The 
choice of a non linear equilibration is made here to high- 
light the the existence of a non linear center manifold 
for the large scale dynamics. In other words, eq.(4) with 
-^^2(^2) = a„ii?2 is supposed to describe the 'normal 
form' dynamics of the large scale magnetic field. Note, 
that our assumption on M2 does not necessarily imply 
a time scale separation between the characteristic time 
scale of i?2 and the magnetic turbulent field. Finally, 
since the system has an inverse cascade of helicity, we 
set i?i = as boundary conditions. 

The free parameters of the model arc the power input 
/o, the magnetic viscosity and the saturation param- 
eters am- Actually, the parameter /o could be eliminated 
by a suitable rescaling of the velocity field. We shall keep 
it fixed to /o = 1 — i. In figure ([T]) we show the amplitude 
of (1-821) and the magnetic energy Eb = (= E„|i?np) as 
a function of Vm for v = le — 7, where the symbol (..) 
stands for time average . For very large v^m the magnetic 
field does not grow. Then, for Vm greater than some crit- 
ical value, {B2) as well as Eb increases for decreasing 
Vm- Eventually, (I-B2I) saturates at a given value while 
Eb still increases, showing that for Um small enough a 
fully developed spectrum of Bn is achieved. This type of 
behavior is i n ag reement with previous studies of Taylor- 
Green fiows ESEI, 32*2 flows in a sphere [111 or MHD 
shell models [23| . In the top insert of the same figure we 
show the magnetic and energy spectrum for = 10^"^. 
Finally, in the lower insert we plot the magnetic dissipa- 
tion = VmX'nk'^{\Bn\'^) and the large scale dissipation 
due to A/„. Note, that at the dynamo threshold, we ob- 
serve a sudden bump in the magnetic dissipation which 
decreases for decreasing Vn. At relatively small Vm, the 
magnetic dissipation becomes constant and quite close to 
the large scale dissipation. 

We can reasonably predict the behaviour of (|i32p) as 
function of Vm by the following argument. The onset of 
dynamo implies that there exists a net flux of energy from 
the velocity field to the magnetic field. At the largest 
scale, the magnetic field B2 is forced by the velocity field 
due to the terms F2{u, B), sec cq.(4). The quantity A = 
TZ[F2(u^ B)B2] is the energy pumping due to the velocity 
field which is independent on B2 and a™. Thus, from 
eq.(4) we can obtain: 



ld\B2\' 
2 dt 



= A- am\B2WBl. - Bi) - ,yrakl\B2\^ (5) 



where B2r and B2i are the real and imaginary part of 
B2. For large Vm, the amplitude of B2 is small and the 
symmetry breaking term proportional to Um is negligible. 
Under this condition, and with the boundary condition 
constrains, we expect from ([5]) or ([7]) that the behavior of 
B2 is periodic, as it has been observed in the numerical 
simulations. On the other hand for relatively small i>rm 
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FIG. 2: Time behavior of B2r for three different values of i/m 
(displayed on the left side) and constant i^. The blue segment 
in the upper panel shows lOOtd, where td is the dissipative 
time scale computed as td ~ 1./ {k^v-m}- One time unit in 
the figure corresponds to the large scale eddy turnover time 



FIG. 3: Average persistence time r as a function of the mag- 
netic viscosity Vm for am = 0.1 and u = 10~^. The green 
line corresponds to the fit given by equation ((9)1 . In the insert 
we plot l./log{r) versus Um to highlight the linear behavior 
predicted by 



the non linear equilibration breaks the U(l) symmetry 
and i?2i becomes rather small and statistically stationary 
solutions can be observed with = A/a„i. Comput- 
ing A from the numerical simulations, we can use ([5]) to 
predict how (|i32p) depends on ;/,„. The results is shown 
in figure [1] by the blue line with rather good agreement. 

We are interested to study the behavior of the mag- 
netic reversal, if any, as a function of and in partic- 
ular in the region where \B2\ saturates, i.e. it becomes 
independent of ■ In figure [21 we show three different 
time series of the i?2r = Re{B2) as a function of time for 
three different, relatively large, values of the magnetic 
diffusivity. The figure highlights the two major informa- 
tions discussed in this letter, namely the obcrvation of 
reversals between the two possible large scale equilibria 
and the dramatic increase of the time delay between re- 
versals for increasing Vm values. Note that this long time 
scale, as observed in the upper panel of figure [21 is much 
longer than the characteristic time scale of B2 near one 
of the two equilibrium states. The system spontaneously 
develops a significant time scale separation, for which 
given polarity is maintained for times much longer than 
the magnetic diffusion time. In figured we show the av- 
erage reversal time as a function of Vm- More precisely, 
let us define t„ the times at which -82 (^n) = and B2 
has opposite sign before and after t„- Then the reversal 
(or persisntance) is defined as Tn = tn — in-i, while the 
average reversal time r is defined as the average of t„ . 

Figure [3] clearly shows that for large i^m , becomes 
extremely large (note that the figure is in log-log scale). 
Thus, even if neither (IB2P) nor Cd depend on Vm, the ef- 
fect of magnetic diffusivity is crucial for determining the 
average time reversal. In order to develop a theoretical 
framework aimed at understanding the result shown in 



figure [21 we assume, in the region where (I-B2P) is inde- 
pendent on that B2i ~ and that the term F2{u, B) 
can be divided into an average forcing term proportional 
to i?2r and a fluctuating part: 

F2{u,B)^ I3B2 + ^' (6) 

where (3 depends on /o and 4>' is supposed to be uncor- 
related with the dynamics of B2 , i.e. ([^'-Bj]) = 0. Note 
that in the context of the mean-field approach to MHD, 
the first term I3B2 would correspond to an 'alpha-effect'. 
Using ([6]), we can rewrite the equations for B2 as follows: 

^^f3B2- arnBl + 4>' . (7) 

where we neglect the dissipative term since /? S> Vmk2 
in the region of interest. Eq.(l7|) must be considered an 
effective equation describing the dynamics of the mag- 
netic field B2 and its reversals, and the fluctuations (j)' 
incorporates the turbulent fluctuations from the velocity 
and magnetic field turbulent cascades. It is the effect 
of (j)' which makes the system 'jump' between the two 
statistically stationary states. Using (O we can obtain 
13 = ^ Attm while the two statistical stationary states can 
be estimated as ±i?Oi B^ ~ (3 /am- The effective equa- 
tion ([7]) is a stochastically differential equation and, using 
large deviation theory, we can predict r to be 







^ ^"P (^) 







(8) 



where a is the variance of the noise (/)' acting on the 
system. Let us notice that A and a must have the 
same dimension, namely [Bj^/time. Thus, we write a 
as a = Af where / is a function of the relevant di- 
mcnsionlcss variables. In our problem the dimension- 
less numbers expected to play a role for the dynamical 
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FIG. 4: Time behavior of B2r for two difFerent values of i/m 
(displayed on the right) obtained by using a linear large scale 
equilibration —7^2 (7 = 0.13) and by imposing B2i — 0. 
Note that although large scale equilibration is achieved by a 
linear damping on the magnetic field, B2r shows quite well 
defined statistical equilibria due to the symmetry breaking 
constrain B2i = 0. 



not able to fit the entire range of Vm shown in figure [3l 

We finally comment on the choice of a non linear term 
in equation (4) . Actually, we can avoid non linear equili- 
bration to obtain the same (qualitatively) results. In fig- 
ure [4] we show two cases obtained with M2{B2) = — 7-B2 
with the constrains i?2i = and 7 = 0.13. The equi- 
libration mechanism is therefore linear while the sym- 
metry breaking is obtained by the constrain B2i = 0. 
Thus the two requirements, large scale dissipation and 
symmetry breaking, are satisfied. Figure [4] shows that 
statistical equilibria can be observed independent of non 
linear mechanism. Moreover, by changing the magnetic 
diffusivity, we can still observe a rather large difference 
in the average reversal time. We argue that this effect is 
independent on the particular choice of the equilibration 
mechanism since it is dictated by dimensional analysis 
and large deviation theory. 
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behavior of the magnetic field are: the Reynolds num- 
ber Rv, the magnetic Reynolds number Rm (or equiva- 
lently the magnetic Prandtl number Pm) and the quan- 
tity Rm = y/ Aum / {'^m^2) which is an effective Reynolds 
number, corresponding to the efficiency of energy trans- 
fers from the velocity field to the magnetic field at large 
scale. Given the fact that we operate at constant power 
input and Ry = const, we expect / to be a function of 
(i?m, Rm) only and we also expect the effective magnetic 
Reynolds number to be proportional to the integral one 
{Rm oc Rm). We then show below that a very good de- 
scription of our numerical results is obtained using the 
lowest order approximation f{Rm,RM) = R\[ ~ Rm, 
where R\,j is a critical magnetic Reynolds number below 
which reversals are not be observed. This choice leads to 
a — A{i'^ — Vm)/uL, and finally to 

r^e^^(-^ ) , (9) 

where C is a constant independent of i^m- This func- 
tional form is displayed in figure [Sj it agrees remarkably 
with the observed numerical values of r for a rather large 
range. In the insert of figure[3]we show 1 /log{T) as a func- 
tion of Vm to highlight the linear behavior predicted by 
aq. ([9]) . The physical statement represented by ^ is that 
the average reversal time should show a critical slowing 
down for relatively large Vm- In other words, we expect 
that fluctuations around the statistical equilibria increase 
as Rm increases. The increase of fluctuations may not be 
monotonic for very large Rm, which explains why we are 



[1] Magnetohydrodynamics and the Earth's Core: Selected 
Works of Paul Roberts, A. M. Soward Ed., CRC Press 
(2003) 

[2] Merrill, R. T., M. W. McElhinny, and P. L. McFadden, 
The Magnetic Field of the Earth, Paleomagnetism, the 
Core and the Deep Mantle, Academic, London, (1996) 

[3] M. Bourgoin, R. Volk, N. Plihon, P. Augier, J.-F. Pinton, 
Aew J. Phys. 8, 329, (2006) 

[4] R. Monchaux et al., Phys. Rev. Lett. 98 044502, (2007) 

[5] M. Berhanu et al., 77, 59007 (2007) 

[6] D. Sweet D et al., Phys. Rev. E 63(6), 066211 (2001) 

[7] Hoyng P, Ossendrijver MAJH, Schmitt D, GAFD, 94(3- 
4), 263-314 (2001) 

[8] A. Gailitis et al., Phys. Rev. Lett. 86, 3024-3027 (2001) 

[9] R. Stieglitz and U. Miiller,Phys. Fluids 13, 561-564 
(2001) 

[10] R. Volk et al., Phys. Rev. E, 73, 046310 (2006) 
[11] Spence et al. Phys. Rev. Lett. 98, 164503 (2007) 
[12] Denisov et al., JETP Letters 88(3), 192 (2008) 
[13] R. Benzi, Phys. Rev. Lett., 95, 024502 (2005) 
[14] L. Biferale, Annu. Rev. Fluid Mech. 35, 441 (2003) 
[15] Frick P and Sokoloff D., Phys. Rev. E57, 4155(1998) 
[16] Stepanov R and Plunian P, J. Turb. 7, 39 (2006) 
[17] Plunian et al., A^'ew Journal of Physics 9, 294 (2007) 
[18] Ponty Y, Politano H, Pinton JF, Phys. Rev. Lett. 92(14), 
144503 (2004) 

[19] Baerenzung J, et al., Phys. Rev. Lett. 78(2), 026310 
(2008) 

[20] Y. Ponty et al., Phys. Rev. Lett. 94, 164502 (2005) 
[21] J-P. Laval et al., Phys. Rev. Lett. 96, 204503 (2006) 
[22] A. Bayliss et al., Phys. Rev. E 75(2), 026303 (2007) 
[23] P. Frick, R. Stepanov, D. Sokoloff, Phys. Rev. E 74, 
066310 (2006) 



